Theory of carrier concentration-dependent electronic behavior in layered cobaltates 
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A natural explanation for the carrier concentration-dependent electronic behavior in the layered 
cobaltates emerges within correlated-electron Hamiltonians with finite on-site and significant nearest 
neighbor hole-hole Coulomb repulsions. The nearest neighbor repulsion decreases hole double- 
occupancy below hole density |, but increases the same at higher hole densities. Our conclusion 
is valid for both single-band and three-band extended Hubbard Hamiltonians, and sheds light on 
concentration-dependent e' g hole occupancy within the latter. 
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Layered cobaltates - anhydrous Na x Co02, LLjCoC^ 
and the misfit cobaltates [Bi 2 A 2 04] [Co02] m , where A = 
Ba, Sr or Ca and m is incommensurate - have attracted 
wide attention for their unconventional metallicity and 
tunability of the carrier concentration. Na x Co02 con- 
sists of edge-sharing Co06 octahedra, with C0O2 layers 
separated by Na layers. The Co ions with average charge 
(4-x) + form a triangular lattice. Both experiments [1| 
and theory Q indicate large crystal-field splitting and 
therefore low-spin states for the Co-ions. Trigonal distor- 
tion splits the tig d-orbitals into two low-lying e' g orbitals 
and a higher a\ g orbital. Charge carriers are S = \ holes 
on the Co 4+ ion sites Q. The hole density p = 1 — x in 
anhydrous NazCo02 and Li :c Co02. Angle-resolved pho- 
toemission from Naa;Co02 indicate that carriers occupy 
the a\g orbitals only although Conipton scattering 

finds small x-dependent e' g contribution [6(. 

The electronic and magnetic behavior of these 
materials exhibit a peculiar carrier concentration- 
dependence. The temperature-dependent magnetic sus- 
ceptibility x(T) in Na x Co02 was early on character- 
ized as "Pauli paramagnetic" for a; < 0.5 and "Curie- 
Weiss" for x > 0.5 0j- The Curie- Weiss behavior reflects 
strong Coulomb repulsion between the holes [H, H, 0, 0, Hi ■ 
Strong correlation at large x is supported by observa- 
tions of charge-ordering (CO) 0, [l(|, Na-ion ordering 
1 1] , spin-density wave and intralayer ferromagnetic cor- 
relations [fjj]- Qualitatively different behavior for the 
small x region is also agreed upon. Recent experimen- 
tal work suggest that (a) x{T) here is weakly antiferro- 
magnetic rather than Pauli paramagnetic, and (b) the 
cross-over between strong and weak correlations occurs 
at x ~ 0.63 - 0.65 rather than x = 0.5 [lB. CO and 
Na-ion ordering are both absent for small x |10j. 

Understanding the x-dependence of the electronic and 
magnetic behavior of Naa;Co02 continues to be a theo- 
retical challenge. It has been suggested that while the 
ai g -only description is valid for large x, holes occupy 
both dig and e' g orbitals at small x [l3j, and the on- 
site Hubbard correlation U is x-dependent. Quantum 




chemical configuration interaction calculations 
larger a\ g — e' g separation than LDA calculations 
(ii) many-body approaches that take Coulomb hole-hole 
repulsion into account 15- lj} do not find the e' g pock- 
ets on the Fermi surface predicted within LDA calcula- 
tions, x-dcpcndence has also been ascribed to differences 
in the potential due to Na layers ll|, I 
tally, Lia;Co02 [IH and Bi "misfit" cobaltates 



Experimen- 
|20| do not 

exhibit ion ordering but nevertheless exhibit very simi- 
lar carrier concentration-dependence 0, [2(| , suggesting 
that this dependence is intrinsic to the C0O2 layers, with 
the interlayer potential playing a secondary role. 

In the present Letter we show that a simple and nat- 
ural explanation of the carrier concentration-dependence 
emerges within ai g -only as well as multiband extended 
Hubbard models. As we will be interested in all three 
families of layered cobaltates including the misfits, hence- 
forth our discussions involve p which is well-defined for 
all systems, instead of x. Following [li} we write the 
three-band Hamiltonian as 

H = - ^ tctpctaaCj?" + A (" le s ~ "'"J 

+ g X! U aP n ia<r n iPa' + V ^ n ianjP- (1) 

iaflaa' (ij)afi 

Here a and f3 refer to the a\ g and e' g orbitals, c\ aa creates 
a hole of spin a on orbital a on site i, ni aa = c\ aa Ci aa and 
n-ia = J2a nia < 7 - 1S t ne nearest neighbor (n.n.) hop- 
ping, A the trigonal splitting, U = Ug'-° and U' = U°%' 
are the onsite intra- and inter-orbital Coulomb interac- 
tions, and V is the n.n. Coulomb interaction. As in 
[ill, we have ignored the Hund's rule coupling based on 
the very small hole occupation of e' g orbitals (see below). 
As both pho toemission experiments and many-body 
theories [HI TBI 17 1 find negligible role of e' g orbitals, we 
discuss the one-band limit of Eq. [T] first. We show that 
the V term is essential within the one-band model for un- 
derstanding the p dependence of the susceptibility. We 
then show that the same effect not only persists in the 



2 



full three-band model, but influences the p-dependence 
of the e' g orbital occupation as well. 

Single-band limit. Terms containing U' and A are ir- 
relevant and U, V, and t aa refer to a\ g orbitals only. We 
write t aa = t and express U and V in units of t. For 
hole carriers < > j|. Existing oj -only theories largely 
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assume U » t [2j or V = 
few studies that have investigated the effects of finite V 
on triangular lattices are either for particular p = 0.5 
or | or use approximate approaches [25[ that 
do not capture the complex p-dependence that is seen 
experimentally and that we find in our exact solutions. 
We consider here realistic finite U and V and investigate 
all p. We do not assume that U and V are p-dependent. 
Rather, we show that p-dependent correlations emerge as 
solutions to Eq. [H 

Two different observations give the appropriate param- 
eter range, (i) At p = 1 Eq.[T]can be replaced by a V = 
Hubbard Hamiltonian with an effective on-site repulsion 
U eS = U - V. Within the V = Hubbard model for the 
triangular lattice, transition to the Mott-Hubbard insula- 
tor occurs for U c s > U c , where U c ~ 5 — 10. [26|. Exper- 
imentally, p = 1 C0O9 is a poor metal close to the Mott- 
Hubbard transition [8|, indicating that U — V < 5 — 10 
for cobaltates. (ii) For V > ^U, p = | would be charge- 
ordered with all sites doubly occupied (Co 5+ ) and vacant 
(Co 3+ ), as shown in Fig. 1(d). The absence of such CO 
indicates V < |Z7. Taking (i) and (ii) together, we con- 
clude that the likely parameter regime is 6 < U < 14, 
1 < V < 4. Our estimate of V/U is close to that of 
Choy et al. [23|. Our estimate of U/t is slightly smaller 

30. 

We now argue that for realistic U and V there occur 
three distinct hole density regions, (i) p < |, where cor- 
relation effects are strongest; (ii) intermediate | < p < |, 
where correlation effects become weaker with increasing 
density, and can be quite weak at the highest p; and (iii) 
p > I , where correlation effects increase again slowly. We 
classify configurations by the number of double occupan- 
cies Nd- Fig. Q] (a) shows the \/3 x V3 Nd — charge- 
ordered configuration that should dominate the ground 
state of p = |. For fixed U, nonzero V creates an energy 
barrier to holes approaching each other. The effective en- 
ergy cost of creating a double occupancy is thus greater 
than U for p < i, which should exhibit strongly corre- 
lated behavior. This situation changes as p increases, as 
seen in Fig. Q] (b), where we have added a single hole to 
the charge-ordered configuration of p — |. The particu- 
lar hop indicated in the figure that creates a double occu- 
pancy costs only U — 3V. There are only three of these, 
and they increase (Nd) very slightly. With further in- 
crease in p, the number of these low energy hops increases 
rapidly, increasing (Nd) ■ In Fig. [ljc) we show the two ex- 
treme charge-ordered configurations for p — ^, one with 




FIG. 1: (Color online) Dominant ground state configurations 
for (a) p = I and (b) one hole added to p = |. Circles labeled 
'h' ('2h') are singly (doubly) occupied sites, with vacancies oc- 
cupying the vertices of the triangular lattice. The n.n. hop 
indicated by the arrow in (b) costs U — 3V. Dominant con- 
figurations for p = |, with (c) Nd = 0, and (d) Nd = iN. 



thus strong mixing of Nd — and Nd > configurations 
for p close to |, even for V < |{7. Adding double occu- 
pancies to the configuration in Fig. [IJd) is prohibitively 
expensive in energy, which implies that for p > | the 
competition of Nd — is no longer with Nd = N max but 
still with Nd — ^N. We expect correlations to slowly 
increase again in this region. 

We have performed exact numerical calculations to 
confirm the above conjectures. As a measure of correla- 
tions we have chosen the normalized probability of double 
occupancy g(p) in the ground state, 
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N d = (Fig. QIc)), the other with N d = N m( 
The configurations are degenerate at U = 3V. 
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g(p)=l and for U = and U — > 00, respectively, for all 
p, and has intermediate values in between. g(p) is thus 
a measure of U e s(p): small g(p) implies enhanced Curie- 
Weiss type x(T) while moderate to large g(p) implies 
weak antiferromagnetic spin-spin correlations [27]. Our 
proposed mechanism suggests that g(p) is small (large) 
for small (large) hole density, provided V is significant. 

We have calculated g for the six triangular lattice clus- 
ters in Fig. using periodic boundary condition. For 
clusters (a)-(c) the number of holes covers the com- 
plete range p < 1. Computer memory constraints restrict 
us to p < 0.88 for cluster (d) and p < 0.6 for clusters (e) 
and (f), respectively. Our calculations are for all realistic 
U and V < |Z7. As the results are qualitatively the same 
in all cases, we report our results for U — 10 only. 

In Fig. [3] we show our results for clusters (a)-(d) for 
U = 10 and V = 0, 2 and 3. Our data points include both 
even and odd Ah, and except for Ah > 14 in Fig.EJd) we 
have determined the total spin 5* in the ground state in 
each case. With few exceptions, S = S m - in , with = 
(i) for even (odd) Ah. S > S m i n is a finite-size effect, 



FIG. 2: (Color online) Clusters investigated numerically: (a) 
N = 12; (b) and (c) N = 16; (d) N = 18; (e) and (f) N = 20 



FIG. 4: (Color online) Same as Fig.[3jfor N = 20, for U = 10 
and V = (circles), 2 (diamonds) and 3 (triangles); (a) and 
(b) correspond to lattices Fig. [2je) and [2jf ) , respectively. 
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FIG. 3: (Color online) Normalized probability of double oc- 
cupancy of sites by holes versus hole density for U=10, V~0 
(circles), 2 (diamonds) and 3 (triangles) on clusters with (a) 
N = 12, (b) and (c) N = 16, corresponding to Fig. HJb) 
and (c), respectively, and (d) N = 18. Points within squares 
indicate S > Smin (see text). 



as for different clusters this occurs at different densities. 
The g-values for the S > S m i n points were calculated 
correctly in accordance with Eq. [2j The dips in g for 
higher S are expected. In every case we have included 
dashed straight lines connecting the neighboring points 
on both sides. The g for S = S m - ln at these points is likely 
bounded by the computed points and the dashed lines. 

In all four cases in Fig. [3J g(p) is nearly independent 
of p for V = 0, but exhibits the p-dependence predicted 
for V 0. The p-dependence is weakest for N = 12. 
In all other cases there occur distinct strongly correlated 
low density region (p < 0.4), where V suppresses g and 
relatively weakly correlated intermediate density region 
(0.4 < p < 0.8), where V enhances g. The predicted 
decrease in g{p) for larger p is also visible in Fig. EJb) 
and (c). 

Fig. [4] shows plots of g(p) for the 20-site clusters of 
Fig. fJl^e) and (f). As in Fig. [3] we have retained the 
points with S > S mm . Distinct density regions (i) and 



(ii), with opposite effects of V are again clearly visible. 
The boundary between strongly and weakly correlated 
regions is p « 0.30, in agreement with recent experiments 
10]. Calculated charge-charge correlations (niUj) (not 
shown) indicate that while near p — \ (N^ = 6 and 7) 
there is tendency to CO there is no such tendency at 
p = 0.5 (iVh = 10). The qualitative agreement between 
our results for six different clusters, along with steeper p- 
dependence with increasing N, strongly suggest that our 
results will persist in the thermodynamic limit. 

Three-band model, p-dependent g is a consequence of 
the competition between U and V at large p 27J and 
is unrelated to dimensionality or frustration. Since ex- 
act calculations are not possible for the clusters of Fig. [5] 
within the three-band model, we have performed three- 
band calculations for a one-dimensional (ID) periodic 
cluster with eight sites, each with one "ai g " and two 
ll e' g " (hereafter a, ei and e%, respectively) orbitals. The 
only differences between ID and the triangular lattice 
are, (i) the Wigner crystal occurs at p = i in ID instead 
of p = i, and (ii) the maximum in g(p) is expected near 
p = | [27] rather than ~. We retain the same U and V as 
for the single-band model. We have taken t aa = t, and 
inter-orbital hopping = 0.1 — 0.3t (a ^ f3), A = 3|f 
and U' 
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0.6U h 



As is common for ID rings, 
we use periodic (anti-periodic) boundary conditions for 
N h = An + 2 (An) 

In Fig. [Sfa), we have plotted g a (p), the normalized 
probability of double occupancy of the a-orbitals by 
holes, within the three-band model (g e (p) varies by less 
than 15% over the entire range of p). The carrier density 
p here is the ratio of the total number of holes and the 
number of a-orbitals, in agreement with the definition of 
p in Naa;Co02. The g a (p)-behavior is nearly identical to 
that of g(p) in Figs. [3] and |H Interestingly, g a (p) be- 
havior is the same for small and large t e , a , in spite of 
moderate hole population n(e) in the e orbitals in the lat- 
ter case. Calculations for smaller A (not shown) indicate 
similar weak dependence of g a (p) on A. 

In Fig. [Sfb) we plot r(e), the fraction of holes that 
occupy the e orbitals, assuming t a p = 0.3£, for (a) non- 
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Conversely, the absence of Na-ion ordering for weakly 
correlated x < 0.5 in Na^CoC^ further suggests that the 
CO and Na-ion ordering are synergistic effects. 

This work was supported by the Department of Energy 
grant DE-FG02-06ER46315. 
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FIG. 5: (Color online) (a) Normalized probability of double 
occupancy of the a-orbitals by holes in the ID 3-band model 
for N = 8, A = 3 (b) The fraction of holes occupying e- 
orbitals. In both, unfilled (filled) symbols are for inter-orbital 
hopping O.lt (0.3i). Diamonds (circles) are for V = (V = 3). 
In (b), the dashed line shows the noninteracting (U — U = 
V = 0) result. For all other points U = 10 and U' = 6. 



interacting, (b) U, U' ^ but V — 0, and (c) V > cases. 
For t a p — O.lt, r(e) < 0.01 and is negligible. Comparing 
the noninteracting and the V = plots, it is clear that 
nonzero U and U' decrease the e' g occupation, the reason 
for which can only be correlation-induced band narrow- 
ing [Til, 17 1 . It is, however, the V > plot that is far 
more interesting: r(e) now shows a peak at the same p 
where g a {p) has a maximum. This shows that correla- 
tion effects due to V play a much more complex role in 
the multiband picture: in the large-p region V reduces 
U e fi (Fig. 03a)) for the same reason as in the single-band 
picture. This reduces the extent of band-narrowing at 
precisely these p, and leads to an increase in r(e) even 
for p '-independent A. Existing discussions of e' g ho. 
cupancy have largely focused on the Fermi surface 



Our work shows that there can occur weak e' g hole occu- 
pation even if only the ai 9 -orbitals form the Fermi sur- 
face. Most interestingly, the e' g occupation estimated by 
Compton scattering experiments shows a p-dependence 
(see Table I in 6]) that is very similar to our results for 

V > in Fig.[5{b), viz., increasing e^-occupancy with in- 
creasing p. This result simply cannot be explained with 

V = 0, as seen in Fig. (5fb). Our calculations show that 
p-dependent x(T) and e' g occupation are manifestations 
of the same many-body effect. 

Summary. Strongly correlated behavior for small hole 
densities and relatively weakly correlated behavior for 
larger hole densities are both expected for nonzero n.n. 
Coulomb interaction. To the best of our knowledge there 
exists no other satisfactory theoretical explanation for 
the observed weakly correlated behavior nearer to the 
Mott-Hubbard semiconducting hole density and strongly 
correlated behavior farther away from this limit. The 
strong tendency to CO at p exactly | and the absence of 
this tendency at p = | arc both understood. The poten- 
tial due to Na-ions, ignored in our work, will strengthen 
the CO even for incommensurate fillings with p < g 11|. 
We do not find CO at p = 0.5, although it is moder- 
ately correlated. Observed CO here 0] is likely driven 
by the cooperative effects of V and the Na-ion potential. 
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